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Abstract 



> 

. The operator mixing matrix for the dimension six operators in the heavy 

^2 , quark effective theory Lagrangian is computed at one loop. The results are 

\ shown to be consistent with constraints from the equations of motion, and 

' from reparametrization invariance. 

o\ ■ 
iu: 

^ ' Heavy quark effective theory (HQET) |[l|,^ is a useful tool for studying the physics of 
hadrons containing a single heavy quark. The HQET Lagrangian has an expansion in 
^ . powers of derivatives divided by the heavy quark mass m, which translates into an expan- 
sion of hadronic quantities in powers of Aqcd/?^, where Aqcd is the non-perturbative scale 
parameter of the strong interactions. 

The HQET Lagrangian can be computed by matching with the full QCD Lagrangian 
at a scale fi ^ m. This has been done at one loop to order for the two-Fermion 

terms in the Lagrangian and at one loop to order 1/m^ for the two-Fermion terms [Q. 

The renormalization group running of the dimension five (1/m) operators in the HQET 
Lagrangian has been computed ^j^- There are several computations of the running of 
the dimension six (1/m^) operators [^-10] in the literature, but the various papers disagree 



with each other. In Refs. ||7|-[lO| the authors do not take into account the effect of a four 
fermion operator which is present in the dimension six operator basis and is related to the 
Darwin term by the equations of motion. A complete calculation including all dimension 
six operators was given in Ref. 0, but we disagree with this calculation in a few entries of 
the anomalous dimension matrix. 

In this paper we compute the running of the dimension six operators of the HQET 
Lagrangian. The coefficients are computed using one-loop running and tree-level matching, 
which makes the calculation particularly simple. We will use the notation of Ref. 0, to 
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make it easier to compare with previous results. The computations are done in a theory 
with one heavy quark with velocity f , and flavors of massless quarks, g^. The covariant 
derivative is chosen to be + igA^T^. With this convention, the HQET Lagrangian to 
order 1/m^ is 

C = -^G'^'^Gt + y qiilpq + Hiv ■ D)h - —JiD^h - ^gha^^C^h 
4 ^^ ^ 2m 4m 

+cdOd + csOs + Yl CiOi, (1) 
where the Darwin and spin-orbit operators are defined as 

Od = -^h {D^Cn v,h, Os = i-£^ha^,{D'^, G^Vp, (2) 

respectively. The remaining operators Oj are four-Fermion operators involving two heavy 
and two light-quark fields. 



2 2 

i i 
2 2 

- ^^YJ'^ ^4qi. of = /^E ^7^^75/1 qa^l^qi. 

i i 

four-Fermion operators involving four light quark fields, 

2 2 

^3 = ^E^i7''?i ?j7M?j, Of = ^Y^i^'^^^qiqjli.l^qj, 

the light quark penguin operator 



(3) 



(4) 



= ^,I:qil^D,G^'^q^, (5) 



and three-gluon operators 



Of = ^.gfABcGiG^'^^G^'^-, 01 = ^^D^ G^^^. (6) 

The identity 

0^j2D^Gi^ G-^ + 2 g Jabc G^G^^^G^, + G% G^^^ (7) 

has been used to eliminate G^^, G^^^ from the effective Lagrangian. There are also 

operators such as h {iv ■ D)^ h which vanish by the heavy quark equations of motion, and 
can be eliminated by a field redefinition of the Lagrangian. Operators O^^ involving four 
heavy quark fields do not contribute in the single heavy quark sector, and will be omitted 
from the Lagrangian in this analysis. 
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The equation of motion for the gluon field, 



D^G^"'' = gvV^'K + 9J2 (l^T^Yq^, (8) 

i 

can be used to further simphfy the effective Lagrangian. It allows one to rewrite Of 
in terms of four-Fermion operators. It also gives the relation 

Od = Of' + -^hT^h hT^h. (9) 

The second term on the right hand side can be omitted in the single heavy quark sector, so 
the equation of motion reduces to 

Od = Of. (10) 

This relation can be used to further simplify the effective Lagrangian by eliminating OiK 
For the moment, we retain O^^ to make it easier to compare the results with Ref. 0. 

The renormalization group scaling of the terms in the HQET Lagrangian involves 

the operators described above, as well as the time ordered product of two 1/m operators. 
The time ordered products will be denoted by 

Okk='-j d'x T [Ok (x) Ok (0)] , 
Okm = tJ d^x T [Ok (x) 0„, (0)] , 

Omm='-J d\ T [0„ {X) Om (0)] , (11) 

where the kinetic and magnetic moment operators are 

Ok = -^hD^h, Om = ■^gha.^G^^'-h. (12) 
2m 4m 

The coefficient of Ok, Ck, is fixed to unity by reparametrization invariance [^. The coeffi- 
cient of the magnetic moment operator satisfies the renormalization group equation MM 



where Ca = 3 is the Casimir of the adjoint representation. We will later also use Cp = 4/3, 
the Casimir of the fundamental representation (not to be confused with cp), Tp = 1/2, the 
index of the fundamental representation, N = 3, the number of colors, and n/, the number 
of light quark flavors. Note that the coefficient of Om in the Lagrangian is —cp- 
The running of the dimension six operators is 
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The zero entries follow because no one-loop Feynman graph exists for that term in the 
mixing matrix. O^^ represents all possible operators involving four heavy quark fields. 
The tree-level matching conditions at the scale \i = m are cd = = 1, and that c^^, 



;f and cf are all zero. The operator mixing equations Eq. (14) then imply that cf and 



stay zero on scaling from m to /x. The values of cb and cs can be obtained by solving the 
renormalization group equations in the 1 — 2 sector. The equation of motion Eq. (|l^) allows 
one to eliminate Of and further simplify the calculation. The only O^' operator that can 
mix with Oy) or Os is Of\ there are no penguin diagrams from the other O, 
The of operators do not mix among themselves, so that is diagonal |12 



f operators. 
Thus elimi- 



nating of using the equations of motion gives a renormalization group equation for cd and 
Cs only in the {O^, O5, Okk-, 0km-, 0mm} sector. The computation of the coefficients is 
straightforward. The renormalization group equation for the spin-orbit coefficient C5 is 
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;ACACkCF- 



The running of C5 is consistent with the reparametrization constraint Cs 
equation can be solved when combined with Ck = 1 and Eq. (O) for cp 
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l/bo 



(15) 
1. The Cs 

(16) 
(17) 



60 = llOyi/S — ATpnf/S is the first term in the /3-function, and we have used the initial 
condition cs {m) = cp ijn) = 1. 

The renormalization group equation for the Darwin term cd is 
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(18) 
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20 32 Cf 
13 



whose solution is 

using the initial condition cd (m) = 1, and Eq. ([1 

Cs (yU) = 2z'-^ - 1, Cf (/i) = Z~^, CD (/i) 



For QCD, Eqs. f|T 
308 r 



(19) 



T9|) reduce to 
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1 



.-13/2 



(20) 



As an example of numerical values, running the 6-quark terms between rrib and rric gives 

Cf = 0.83, Cs = 0.65, cd = 1.57, (21) 

where we have used nj = 4 and as {rric) jois {rrib) = 1.7, to be compared with the tree-level 
values Cf = Cs = Cd = 1- 

To compute the full effective Lagrangian, including the four-Fermion operators, one 
needs the anomalous dimension matrix Eq. (|T^) . The matrix is computed without using the 
equations of motion to eliminate O^', to make it easier to compare with earlier calculations. 
The entries of the anomalous dimension matrix are listed below. The second form of the 
matrix uses the explicit values of Ca, etc. for QCD. In deriving these equations, we have 
used the identity 

^ {t^, T^}iJX {t^, T^} X = (l - 1/iV') XX + {N- 4/iV) ^jT^^^j xT\ (22) 

which is valid for Fermions in the fundamental representation of SU{N), and holds regardless 
of the 7-matrix structure of the fermion bilinears. 
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12^^ - 2bo = 14 + 4n//3. 
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In writing Eq. ([I^) , penguin diagrams from O^' have been rewritten in terms of O'' and O^^ 
using the gluon equation of motion. The Lagrangian in the single heavy quark sector does 
not depend on Fjs, which have not been computed. 

The anomalous dimension matrix Eq. (|Tip has been computed previously. Fn was com- 
puted in Refs. [[^-[TO 



The submatrix in the 3 — 4 sector was computed in Ref. |]6|,p!3|- |15|] . F, 
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was computed in Ref. • The rest of the matrix was computed in Ref. [§[ . We disagree 

with the previous computations in a few entries of the anomalous dimension matrix. One 
can check the consistency of Eq. (14) with the equation of motion Eq. (|T0|) . One finds, 
neglecting O^'^ operators which were also neglected in Eq. (p!0|). 
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so that the equation of motion is multiplicatively renormalized. This is consistent with the 
result that equations of motion can only mix among themselves under renormalization |]18 . 
One can then use the equations of motion to eliminate O^K The resulting renormalization 
group equations have the same form as Eq. (|T^), with omitted in the second operator 
block. The new entries of the operator mixing matrix will be denoted by T[-. The resulting 
matrix is (note that the zero entries are not the same as Eq. f|]T 
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where the new entries are 
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One can solve the renormalization group equations Eq. 
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to obtain the coefficients in the effective Lagrangian 



with one-loop running, and tree-level matching. The results are Eq. 
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One can solve Eq. (|1^ for cd and cj^' separately. The result for the sum cn + c^^ is the same 
as the value oi cd in Eq. ([I9l), where the equation of motion was used to replace hy Od- 
An independent linear combination of cd and df' that has a relatively simple scaling law is 
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The numerical values for the coefficients with the same choice of parameters as before are 



„hl 
^2 



-0.19, = 0, cf = -0.09, CD = 1.61 and c^^ = -0.04. 
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